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Abstract 

A set of regularly distributed transmission eigenvalues generates a density function. We use such 
a density function inversely determines the form of the indicator function. Using the entire function 
theory, we reduce an uniqueness problem with interior transmission eigenvalues induced by time- 
harmonic Maxwell equation to an uniqueness problem in entire function theory. In such an inverse 
problem, the definite integral of the square root of refraction index is the main parameter. 
MSC:35P25/35R30/34B24/. 
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1 Introduction and the Main Result 



In this paper, we consider the time-harmonic Maxwell equation with non-absorbing refraction index in 
the following setting: 



n(x) = n(r) > 0, r = \x\, when r <E [0, a]; 



0; n(r) = 1, when r > a; n e C 2 (I 



(1.1) 



such that 



V x E t - ikH 1 = 0, 

V x E - ikH = 0, 



V x Hi + ikn(r)Ei = 0, in B; 

V x H + ikE Q = 0, in B; 



(1.2) 



with boundary condition 



v X (Ei - E ), vx(Hi- H ) = 0, on dB, 



(1.3) 

where Eo, Ho is an electromagnetic Herglotz pair, B is an open ball of radius a in M 3 with exterior unit 
normal vector v. We will look for a non-trivial solution to this homogeneous electromagnetic interior 
transmission problem (1.2 1 and (1.3). For each k € C such that (1.2) and (1.3) has a set of non- 
trivial solution is called an interior transmission eigenvalues. We reduce such an electromagnetic interior 
transmission problem to the acoustic interior transmission problem: 



Aw + k 2 n(r)w = 0, 
Av + k 2 v = 



dv 



c)v 
Or 



in B; 
in B; 
on dB, 



(1.4) 



where ii),u6 C 3 (B). To see this, we consider the following quantity 

Ei(x) := V x {xw(x)}; 
ffx(:r) := iV x {E 1 (x)}; 
E (x) := V x {xv(x)}; 
H (x) := iV x {E (x)}, 



(1.5) 
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from which one can obtain a set of solution to the electromagnetic interior transmission problem (1.2) 
and (1.3). We refer the induction to the Colton and Kress [3]. 

We need to consider the solutions w,v to (1.4) that are not spherically symmetric. Therefore, we 
look for non-trivial solutions w, v in the following form: 



v{r,9) = a l j l (jr)P l (cos9); 



w(r, 9) 



-P/(cos( 



(1.6) 
(1.7) 



where Pi is Legendre's polynomial, ji is the spherical Bessel function of degree I, a; and bi are constants 
to be determined and the function yi is a solution of 



(k 2 n(r) 



1(1 + 1) 



)vt = 0; 



]im{^--j l (kr)} = 0. 

r— ¥0 r 



(1.8) 
(1.9) 



for r > and yi is continuous for r > 0. Moreover, as demonstrated in [3], we consider the non- 
spherically symmetric w,v. In such a magnetic problem, we are asked to consider I > 1. Furthermore, 



we see that (1.9) implies 



W(0)=0; wj(0) = 0. 



(1.10) 



We will show there exist a set of k € C with its maximal density and constants ai = a/(fc), bi = bi(k), 
such that (1.6) and (1.7) is a set of non-trivial solution to the interior transmission problem \1A\ . 
Considering (1.5), we see that, for any such value of k, the set of the electric far field patterns is not 
complete in certain functional space. See the discussion in [4]. 

The the interior transmission problem ( 1.4 ) and ( 1.5 1 admits a set of non-trivial solution v, w if there 
exists a set of non-trivial solutions ai , bi to the following homogenous system 



k aiji(ka) = 0; 

a 

ar r 



aikj[(ka) = 0. 



Such a system admits a set of non-trivial solutions a; , bi if and only if the determinant 



di(k) := det \ 

\ dr ^ r 



-ji(ka) 
-kj[(ka) 



(1.11) 
(1.12) 

(1.13) 



di(k) is an entire even function. See [5]. 
In this paper, we consider 



1 = 1; 

di(k) := D(k), k e C; 
y 1 (x, k) := y(x, k), k E C. 



(1.14) 
(1.15) 
(1.16) 



May we ask that if the set of the interior transmission eigenvalues of (1.2), in particular, the set of 
interior transmission eigenvalues of the acoustic system (1.4) or zeros of D(k), can uniquely determine 
the refraction index n(r)l 

Following the local uniqueness results in |12[ 113], we state the uniqueness result in this paper. 



Theorem 1.1 Let the junctional determinant D(z) be defined as in (1.13) and (1.15). Then, the zeros 
of D(z) inside any of the angular wedges 



51 := {k e C| - e < argfc < e}, 

5 2 := {k e C| 7T - e < argfc < ir + e},Ve > 0, 



(1.17) 
(1.18) 



uniquely determine between the spherical refraction indices n(r) if they have the same value at r = 0. 
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2 Preliminaries: Cartwright's Theorem Versus Wilder's Theorem 



Since 



From (1.13), we compute the D(k) as follows. 

D(k) = (-)y(a, k)j[(ka) + 2 ^ a ^j 1 (fca) - \y(a; fc)jx(fca). (2.1) 
a a a z 

. . smt cost ,„ „. 

3i(t) = -p r , (2-2) 

nfn ,-fc, , ,..2cos(fca) ,1 2 . 

D(fc) = (_) p(fllfc ) [ - 1 ^ + (_-_) ffllljfco] 

j/(a, fc) .sin fca cosfca. 1 . sinfca cosfca 
+ a ^Jkaf " ~~fc^ J ~ ^ y(a; ~ ~fc^ J 

,— fc. , , . .cosfca .1 1 . . , . y'(a, fc) .sin fca cosfca. ,„ „. 

= ( — )y(a;fc)[ r — + - smfca+ »i±i . (2.3) 

a (afc) z fca (fca) J a (fca)^ fca 

To study the asymptotics of D{k), we study the asymptotics of the fundamental solutions y(a,k) 
and y'{a,k). For this we use the asymptotic expansion theory for Sturm-Liouville problem provided in 
Erdelyi [8] on page 84 with notations there. 

q(r,k) =n{r)k 2 - \, (2.4) 

which is considered as a power series in -r- In particular, we set 

q (r) = n(r), 9l (r) = 0, q 2 (r) = (2.5) 

Surely we have n(r) > 0, Vr > 0, but 92 (^) is singular at r — 0. Moreover, to fulfill the requirements in 
[8], we use 

AT = 1, S := {fc € C|R{fc[-fl,(r)]*} ^ 0}. (2.6) 

In this paper, we take 

S = {fc e C|0 < argfc < vr} U {fc e C| tt < argfc < 2tt}. 



Accordingly, we have a set of fundamental solutions yi(r), y%{r) to the problem (1.8 1 such that in the 
sectorial region S 

y j (r;k)=Y j (r)[l + 0(±)}; (2.7) 
y' j (r;k)=Y!(r)[l + 0(±)}, (2.8) 
as fc — > 00 in S 1 , uniformly for £ < r < a and for argfc, where 5 is arbitrarily small and 

Y 3 (r) = exp{/? oj fc + 0y + £2^}, (2.9) 



where /3cy, and /?2j satisfy 



(£' ) 2 + n(r) = 0; (2.10) 



Vy + ^ = 0; (2.11) 
2i8oA--^ + W a + ^y=0. (2.12) 
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Following from (2.5), 



Atf(r) 
Pi 3 {r) 
Mr) 



±* / Vlo(p)dp + C j-, 
Jq+ 

--ln[q (r)]+C ir , 



2 
f 2 



dp + a 



0+ 



2.; 



-(fl») a -|8£(p) 



±i/3 2 (r) + C 2j -; 



(2.13) 
(2.14) 
(2.15) 

(2.16) 
(2.17) 



where Ccy, Cy, Caj are constants. Since n(r) > over [0, a] and is in C 2 (R), (3oj( r ), flij(r), fhj{r) are 
all sufficiently smooth and bounded over [0, a]. That is we can set 



/9a(r) 



o+ 



2y / n(p) 



dp. 



In addition, we see 



#i(r) 



4<to(r)' 

l^(r)g (r)-(^) 2 (r) 
4 



Therefore, 



/3 2 (r) 



dp. 



In general, for n > 3, j = ±1, 

M»0 



2^/r^) 

2^(r) 
2 vMp) 



dp. 



(2.18) 

(2.19) 
(2.20) 

(2.21) 
(2.22) 

(2.23) 
(2.24) 



Now we use p.13 ), ( |2.14[ ) and p.17 ). The general solution to ( |1.8 ) is of the form 

y( r ) = [ a - _ e ( i / ''\/9o(P) d ' > ) fc + i ' 3 2(r)i + ^ 1 j-r ^ (p)dp)k-i0 2 (r) i j ^ + (2.25) 

[go (>")]* [(toO")] 3 fc 

We consider a = a{k), /3 = (3(h). To find such a set of constants, we compare ( 2.25[ ) with its asymptotic 
behavior on Oi + K. For this, we combine from [1] p. 263] with an errata in [5, p. 3] that 



y(r) = 



k[n(r)] 1 



(2.26) 



That is a = 47r ; (3 — j^e™ . Since y(r, k) is an entire function in k, (2.25) and (2.26) combine to give 



y(r,k) 



k[n(r)] < 



- cos [k / [n{p)\ 2 dp — 7T + 



/3 2 (r) 



; J [i + o(-)],c\oi+: 



(2.27) 
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uniformly in r, arg/c. Hence, (2.8) gives 



y'(r,k) = {-[n(r)]s 



; ^ r) -r}Bin[* / Hppdp-7C+^][l + 0(j)} 



fc 2 n(r)< 



Mr) 
k 



4k ' 



Mr), 



(r) - * n'(r) cos [k / [n(p)]*dp - tt + ^^][1 + 0{-)] 



(2.28) 



in C \ Oi + : 



//'(..//)= -l+-^]sm[/.- / [n{p))^dp--K+^][l + 0{\)] 



k 



(2.29) 



Using this, we compute the functional determinant D(z) in ( 2.3 ) in C\0i+K. We obtain from ( 2.3 ), (2.27 ) 
and K 28b that 



D(z) 



We set 



Therefore, we have 



cos [z I [n(p)] 2 dp — 7r + ^ 2 ^ ] cos [za] 



^][! + 0(-)] 
z u^z^ z 

-<-,.s[:/ ;i - + Ml ]sinH[ zi + _L] [ i + o(i) 

z a z z a*z J z 

/3 2 (r) 1 
+ sin[z / [n(p)] 2 dp — tt H ]cos[za][-^— 

+ sin[z / [n(p)]4rfp-7r+^]sin[za][^ + ^][l + 0(-)] (2.30) 

z a a z* a°z' i z 



cos z(a — / [n(p)] 2 c?p) + 7r 



Mr) 



}[i + o(- z )} 

Mr), 



+ [-727 + 7473]Nnz(a- / [ n (p)]*dp)+n-^][l + 0(-)] 



sin{z J^[n(p)} Up-TT+ ^} sin[za] ,1., 
+ a 5 z 4 L ^ ^ JJ " 



iV(a):= / [n(p)]*dp. 



(2.31) 
(2.32) 



D(z) 



-r^[cos za - zN (a) + tt - + 0(- 
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-!-[ -^ + ^][sinza-ziV(a)+7r-^][l + 0( 1 )] 
arz a*z J ■ 



sin{zjV(a)-7r+feM}sinH ri , n ,l, 
IsTi i 1 + ^"Ji- 



(2.33) 



Actually, it has been shown in [3] if considered with the errata in [5J p. 3] that 

1 



di(z)\ i=1 = -^{sin{z(iV(a) - a) - tt} + 0( )}, z € 0« + ! 

arz 



, lnz, 

z 



(2.34) 



We see (2.33) matches with (2.34) when limited on Oi + E. Without loss of generality, we discuss the 
zeros of 



z 4 D(z) := D t (z) + D 2 (z) + D 3 (z), 



(2.35) 



where Di{z) is the i-th term in (2.33) with a multiple of z , i = 1,2,3. For this representation form of 
an entire function, we consider the value distribution theory from complex analysis which we refer to the 
Levin's book PHE] and Cartwright's theory [2J[3]. Let us review the following definitions. 
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Definition 2.1 Let f(z) be an entire function. Let Mf[r) :— max| z | =r |/(z)|. An entire function of f(z) 
is said to be a function of finite order if there exists a positive constant k such that the inequality 

M f (r) < e r " (2.36) 

is valid for all sufficiently large values of r. The greatest lower bound of such numbers k is called the 
order of the entire function f(z). By the type a of an entire function f(z) of order p, we mean the 
greatest lower bound of positive number A for which asymptotically we have 



That is 



M f (r)<e Ar ". (2.37) 
ff:=limsup ^/M. (2 .38) 

r— >oo r? 

If < a < oo, then we say f(z) is of normal type or mean type. 
We also see that 

e (°-t)r" < M f{r) < e^+^ r ", (2.39) 

n as 

where the first inequality holds for some sequence going to infinity and the second one holds asymptoti- 
cally. 

Definition 2.2 If an entire function f(z) is of order one and of normal type, then we say it is an entire 
function of exponential type a. 

Definition 2.3 Let peR and p(r) : K + — > M. + . We say p(r) is a proximate order to p if 

lim p(r) = p > 0; lim rp'(r) lnr = 0. (2.40) 

Definition 2.4 Let f(z) be an integral function of finite order in the angle [9\ , 82] ■ We call the following 
quantity as the generalized indicator of the function f(z). 

In I f(re l9 )\ 

^(g) : =limsup U y )l ,e 1 <9<9 2 , (2.41) 
where p(r) is some proximate order. 

The order and the type of an integral function in an angle can be defined similarly. The connection 
between the indicator hf(8) and its type a is specified by the following theorem. 

Lemma 2.5 (Levin [10j,p.72) The maximum value of the indicator hf (9) of the function f{z) on the 
interval a < 9 < j3 is equal to the type a f of this function inside the angle a < argz < j3. 

Lemma 2.6 Let A, B , C and D be real constants. 

h cos {A z +B}(0) = \Asm6\; (2.42) 
h sin{C z+D}(0) = \Csmd\. (2.43) 



Proof We observe that 



I cos Az + B\ 2 = sinh 2 {Ay} + cos 2 { + B}; (2.44) 
I sin Cz + D\ 2 = sinh 2 {C*y} + sin 2 {Cx + D}. (2.45) 



Applying definition (2.411, we prove the lemma. □ 

Lemma 2.7 For any entire function D(z), the following identity holds. 

h z 4D(z){9) = h D ( z) {9). 



G 



Proof z i D{z) and D(z) are both entire function. We see that 

ln|z 4 J(z)| 41n|z|+ln|D(z)| 
« z 4 D ( 2 ) (9) = lim sup = hm sup = h D ^. 



)(<?)• 



(2.46) 



□ 



We mention two more inequalities for indicator functions. 

Lemma 2.8 Le£ /, g &e £wo entire functions. Then, the following two inequalities hold. 

hfg(6) = hf(6) + hg(0), if one limit exists; 
hf+g < max{h f (6),h g (6)}, 



(2.47) 
(2.48) 



'2.35) 



where if the indicator of the two summands are not equal at some 9q, then the equality holds in (2.4-8). 
Proof We can find these in [TOl p. 51]. □ 

Lemma 2.9 If a ^ \J n(p)dp, then the indicator function of D^{z) 

h D3 (9) = (a + N(a))\smO\; 

h Dl {9) = |JV(a)-a||sin0|; 
h D2 {9) = \N(a)-a\\sm6\; 
h 7 A D {0) = (a + N(a))\ sin 9\. 



(2.49) 



// a = J a \J n(p)dp, then 



h z i D (9) = 2a| sin^l 



(2.50) 



(2.51) 



Proof The indicator function of cos{zN(a)} is |JV(o)|| sin0| and the one forsmza is a|sin0|. Apply- 
ing (2~J1^ and (2~P\ ) to D 3 (z), we prove the (2.4Sty - Similarly, we use (2.42\ ) and (2.43\ l to prove the 
ho 1 (9) and hn 2 (9) respectively. Accordingly, we apply (2.48 ) to obtain 



L Dx + D. 



m&x{h Dl (9), h D3 (0)} 

i{\N(a) - a 
\N(a) + a\\ sin( 



max{|7V(a) - a|| sin0|, \N(a) + a|| sin0|} 





Finally, using (2~J^ and (248} again, 



ha(6) — h.D 1 +D 2 +D 3 (9) 
= max{/i Bl+ D 3 (e), h D2 (9)} 



= (o + iV(o))|sin0|. 



(2.52) 
(2.53) 



(2.54) 
(2.55) 



To prove (2.51), we begin with { 2.33]) . 



lim 

r— f oo 



In | sm{zN(a)} sin{za}| 
r 

In \{cosz(N(a) - a) - cosz(N(a) + a)}/2\ 



(2.56) 



Since a — ^n(p)dp, (2.51) follows from (2-42). □ 

Definition 2.10 We say an entire function is of class C if the following two conditions hold: 

f(z) is an entire function of exponential type; 
log + |/WI 



1 + t 2 



-dt < oo. 



(2.57) 
(2.58) 
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Lemma 2.11 The entire function z i D{z) in (2.35) is an entire function of class C. 

Proof The condition ( 2.58| ) is direct from (2.34 1. From (1.13), each term there is an entire function of 
normal type. Hence, z A D(z) is of class C. □ 

We need more tools to understand the distribution of the zeros of z A D{z) and its analytic behavior. 

Definition 2.12 The following quantity is called the length of the indicator diagram of entire function 
/• 

d=h f C-) + h } {-*-). (2.59) 



Lemma 2.13 // (2.51) and (2.58) hold, then f(z) is of completely regular growth and all of the zeros of 
f(z), except possibly a set of zero density, lie inside arbitrarily small angles | argz| < e and | argz— ir\ < e, 
where the density 

Al ;= Km ^ (2.60) 

t— > oo t 

of the set of the zeros within each of these angles is equal to d/2ir, where d is the length of the indicator 
diagram, i = 1,2 denotes the quantity corresponding to each of the two such angles along the real axis 
and Ui(t) is the number of zeros in each angles with norm less than t. 

This theorem is in Levin (TUJ p. 251]. 

Lemma 2.14 The length of the indicator diagram of z 4 D(z) is 2(a + N(a)) . The density in each of the 
two small angles along real axis is (a + N(a))/iT. 



Proof This follows from ( |2.50[ ) and definition ( |2.59[ ). □ 

Such a density theorem is crucial in this paper. We present an alternative justification from the point 
of view of the zeros of exponential sums. Let's discuss Wilder's type of theorem for exponential sums. 
We refer to the work of D.G. Dickson [HI [7] and R.E. Langer [5]. The reference [Hj is a concise review 
and is a detailed lecture. The reference [3] is a step-by-step introduction to such a theory. 

Let 



f{z) = Y,A j z mi [l + e{z)]e^*, 

3=1 



(2.61) 



where n > 1 and Aj and ujj are complex numbers such that Aj ^ and the OJj are distinct; the rrij are 
non-negative integers; the functions e are analytic for \z\ > ro > with lim 2 _ yoo e(z) — 0. When we are 
talking about the zeros of f{z), we are referring to its zeros outside certain open ball around the origin. 
A function of the form 



3=1 



(2.62) 



where Pj{z) are polynomials, is a straightforward example of (2.61). Our z D(z) in (2.35) would yield 
another special case: 

Z i D(z) = \ Cie -iN{a)z-iaz + ^jN^z-iax + c ^ e iaz-iN{a)z + W(o)«+ia»l M + 0(-)], (2.63) 

z 

where C2, C3 asymptotically approach to constants as z — ¥ 00 and cy, C4 are constants. 
An introductory version of Wilder's theorem is of the following form. 

Theorem 2.15 (Dickson) Let 

R(a, s, H) :={z = x + iy e C| |a;| < H, y e [a, a + a]}; (2.64) 

N g (R(ot, s, H)) := { the numbers of g{z) in R(a, s, H)}. (2.65) 

Let g(z) — y^j—i Aje^i* , where z — x + iy, Aj ^ ; u>i < aj 2 < ■ ■ ■ < w n . Then, there exists K > such 
that 
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1. each zero of g is in \x\ < K; 

2. for each pair of reals (a, s) with s > 0, 

\N g (R{a,s,K))-s(u n -LJ l )/(2n)\ <n-l. (2.66) 

We refer these to [BJ. 

To consider our modified determinant z 4 D(z), we need a more sophisticated version of such a theorem. 
We set up the following quantities to the f(z) in (2.61 ): let Q be the broken line given by the uJj given 
in (2.61) with wi, ■ • • ,uj a as its vertices. The indices is labelled counterclockwise. Let Lk be the line 
segment p k ,Uk+i] and 

<f>k ■= arg{u; fc - oJ k +i} 

in [-§ , ^f). Let 

Certain tJ p on L^ are assigned doubly indexed subscripts as follows: let the convex hull of tJfc, aJ^+i and 
t p = UJ p + im p e k in which Zo p on assign subscripts j = 1, • • • , Ofc to Wjy so that = u^, w^o-,. = Wfe+i 
and Tkj are vertices of this convex hull and preceding in a counterclockwise direction from TD k + im k e k 
to oJfe+i + im k+1 e k . For j = 1, • ■ • ,cr fe - 1, 



Lkj '■— [ T kj, T kj- 



1 



m fci - m fc j+i 

Mfej — 7 s — (2-67) 

which is real; n k j is the number of t p on Lj.,-. 

Moreover, for j = !,■■■ ,o k — \ and H > 0, we define 

V fci (JT) := {«| 3(z/e fe ) > 0, |R(z/e fc ) + /i fci log |z|| < ff}; (2.68) 

Tk(9) is defined to be a closed sector with vertex at zero of opening 26 about the outward normal to L k 
through the origin. For the same k and j and each triple of reals (a, s, H), s > and H > 0, the set 

R kj (a,s,H) := {z|3(z/e fe ) +^ fcj -argz € [a,a + s], |^(z/e fe ) + /x fcj - log \z\\ < H}, (2.69) 

where arg z € 0& +tt). It's a curvilinear tubular neighborhood along some direction. We may collect 
a few facts from [51 [7] . 

Lemma 2.16 The following properties hold. 

1. The size of R k j(a, s, H)is approximately a rectangle of dimension s by 2H for large a. 

2. V k j is individually connected and disjoint to each other for each pair of (k,j) when \z\ is large. 

3. The boundaries of V k j are logarithmic to the outward normal to L k to the exterior ofQ. 
4- For each fixed 9, H, the subsets ofV k j(H) are in T k (9) for large z. 

5. For fixed 9 and H , 

R k j(a, s, H) C V kj n T k (9), Vs > 0, for large a. (2.70) 

Let us define 

N f (R k j{a, s, H)) := { the numbers of f(z) in R kj (a, s, H)}, (2.71) 

where the zeros are counted according to their multiplicity. We refer to the page 4 and 10 in [7] for some 
delicate diagrams on all of these quantities. Most importantly wc have the following theorem. 



Theorem 2.17 (Dickson) Let f{z) be given as in (2.61). Then, there exists K > such that (a) all 
but a finite number of zeros of f of modulus greater than r$ are in Uk.jVkj and (b) for each pair of 
positive reals e and so, there exists an c<q — ao(e, sq) such that whenever a > ao and s > sq, 

\N f (R kj (a,a,K))-8\u kj+1 -u kj \/{2ir)\ <n kj -l + e. (2.72) 
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Figure 1: The polygonal Q 

This is stated as in [6]. 

Without loss of generality, we consider a Wilder's type of theorem to (2.63) in the following setting: 
V{z) := {z/i) 4 D(z/i) = c ie - N{a)z - az + c' 2 z 3 e N{a)z - az + 4z 3 e az - N{a)z + c 4 e N{a)z+az (2.73) 

and 

a > N{a). (2.74) 

We set 

loi = uJi = a + N(a), uj 2 = uj 2 = —a — N(a), as points in C. (2-75) 
Therefore, we do the double-indexing for T)(z) as follows. 

Un = Tn = a + N(a), uj 12 — a — N(a), w 13 = —a + N(a), aJ 14 = t 14 = —a — N(a). (2.76) 

We refer to Figure 1 below for a diagram. Furthermore, the line segment 

Li = [wi4,wii], Lu = [wn,Ti 2 ], L12 := [ri2,ri 3 ], L 13 := [ti 3 ,cju]- (2.77) 

These are edges for polygon Q. 

(j>i = arg(wi - cj 2 ) = 0. (2.78) 

Henceforth, we have 

ei = 1. (2.79) 



By (2.73) and the dchnition of t p , 

mu = 0, mi2 = 3, mi3 = 3, mu = 0. (2.80) 

Therefore, 

Tl2 = a- N(a) + 3i; t 13 := N(a) - a + 3i. (2.81) 
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We compute more. From (2.67), we obtain 

-3 



2iV(o) ;Ml2 ~° ; All3_ 2iV(a)' 



Accordingly, we have 



V n ={z\$S(z)>0, 
{z|3(^)>0, 
y 13 = {z|3(z)>0, |K(z) + 



27V(a) 

V 12 = {z|3(z)>0, |R(z)| <iJ}, 

3 

2iV(a) 



log|z|| < 
log|z|| < H}. 



(2.82) 

(2.83) 
(2.84) 
(2.85) 



Combing (2.731 to (2.85) and Theorem 2.17 we prove the following theorem. 
Theorem 2.18 There exist some K > and large a such that 

S ^M, a>N(a); 



N v (R n (a,s,K)) 



N v (R 12 {a,s,K)) 



s%, a<N(a); 
a-N{a)\^ 



L s~, a< N (a). 



(2.86) 
(2.87) 
(2.88) 



For the critical case a — N(a), we have to reset the polygonal Q. We define: th = uj\ = uJn; oJn = 
Tn = 2a, U)\2 = 0, a7i3 = T13 = —2a. Accordi ngly, ran — 0, mi2 = 3, mi3 = 0. Moreover, /in = — 
/i 12 = For a quick picture, we just delete (2.87) for the critical case. 



Corollary 2.19 Let Aj, £ = 1,2, &e defined as in Lemma \2.13\ There exists some ao such that 

A, = N v {U 3 j= iRi 3 (a,s,K)) = whenever a > a . (2.89) 



Proof Since {Vfej} are disjoint to each other for large a, we sum up all three asymptotic behaviors 
from ( |2.86[ ), ( |2.87[ ) and ( |2.88[ ) to conclude that 



a + TV (a) 



(2.90) 



Also, [wi4,wi3], [wi3,wi2], [wi2,wn] have the same outward normal, so i?n(a, s, K), Ri2(a, s, K) and 
i?i3(a, s,if) asymptotically fall in the same wedge T(9) around this outward normal with any fixed 9. 
Let 9 be small. We prove the corollary. □ 



We see that (2.84) and (2.87) is weaker than the transmission eigenvalues distribution theory in [14] . 
However, (2.86) and (2.88) offer some windows outside the real axis. 

3 Proof of Theorem 11.11 

Let D t (z) be the functional determinant corre sponding to refraction index n l (r), i — 1,2. If the zeros of 
D l in either wedge coincide, then Lemma 2.14 tells us that the indicator diagrams from the two refraction 
indices, say, cf, are identical. That is 



d 1 :=h D i(^) + h m (-~)=d 2 . 



Then, from (2.50) and (2.51), 



h D i{9) = h D 2(9), V0. 



In particular, let N l (a), i=l,2, be defined as in (2.32). Lemma 

jV 1 ^) = N 2 {a). 



2.14 



says 



(3.1) 
(3.2) 

(3.3) 
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Let y l (r, z) be the asymptotic solution given as in (2.27) with respect to index n l (r). Now we apply (3.3) 
with ( |2.27 1 to the following entire function of normal type, 



so we have 



F(z) :=y 1 {a,z)-y 2 (a,z),C\Oi + : 



, f rw \n\y 1 (a,re ie )-y 2 (a,re ie )\ „ n , „ 
h F (9) = limsup ly v ' - — - = 0, (9 ^ 0, tt. 



Since hp(0) is continuous in G [0, 2tt], we have 

h F (0) = 0. 



(3.4) 
(3.5) 

(3.6) 



Therefore, using Lemma 2.5 we conclude that F{z) is an exponential function of zero type. 

As inspired by Aktosun, Gintides and Papanicolaou [1 , McLaughlin and Polyakov [12] , McLaughlin 
[13[ H4| . we seek to identify n(r) by inverse Sturm-Liouville theorem. In particular, we apply the results 

in eng. 



Theorem 3.1 (McLaughlin) We consider the following Sturm-Liouville problem 

z" + (k 2 - q)z = 0, 0<x<l; 
z(0) = z(l) = 0, 

where q G L 2 (0, 1). For another boundary condition, 

z(0) = z'(l) = 0. 



(3.7) 
(3.8) 



(3.9) 



Suppose gi, q2 G L 2 (0, 1) and, X n {qi) — Xn(Q2)> the eigenvalues to (3.1) and (3.8), /i n (<Zi) = Mnfe), the 
eigenvalues to (3.7) and (3.9), Vn G N. Then, q\ = q%, a.e. 

We need more lemmas to apply such a theorem. One is the following Phragmen-Lindelof type of theorem 
in Levin [11[ p. 38, Theorem 3]: 

Lemma 3.2 If f(z), z = x + iy, is an analytic function in the half-plane {z\ 3z > 0} such that, for all 

€ > 0, 

M f (r) < e ( CT+£ ) r , 



and | f(x) | < M on the real axis, then 



\f(x + iy)\ < Me ay . 



Lemma 3.3 Letn l {r), i = l,2, be two refraction indices satisfying (1.1) generating same set of interior 
transmission eigenvalues in either Ei or £2 ind n 1 (0) = n 2 (0). Let y l (a,z) be the solution defined by 
index n l {r). Then, y 1 (a, z) = y 2 (a, z) and {y l )'(a, z) — (y 2 )'(a, z) = 0. 



Proof The type of y (a, z) — y 2 (a, z) as an entire function of exponential type is zero by Lemma 
That is a = for previous lemma. Moreover, given iV : (a) = N 2 (a) and (2.26), 



2.5 



F(z) = y 1 (a, z) — y 2 (a, z) — > 0, as z — > Oi ± 00. 



(3.10) 



Using Lemma 3.2 a nd Li ouville's theorem in complex analysis, we conclude that y 1 {a, z) — y 2 (a, z) = M, 
for some constant. ( 3.10 1 again implies that y 1 (a, z) — y 2 (a, z) = 0. 
Similarly, from (2~29), we obtain (^/(o,^) - (y 2 )'{a,z) = 0. □ 



Now we consider the Liouville transformation of (1.8 



z (0 [n-ir^Viir), where £, := J^[n(p)]^dp. 



(3.11) 



In this case, (1.8) becomes 
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where 



P(0 := 



z" + [fc 2 -p(£)]z = 0; 
*(0) =0; z(JV(a)) = y(a). 

n"(r) 5 [n'(r)] 2 



(3.12) 



(3.13) 



4[n(r)] 2 16 [n(r)] 3 r 2 n{r) 
The zeros of y(a, A:) exactly corresponds to the eigenvalues of the Sturm-Liouville problem 

f y" + (fc 2 n(r) - - 0, < r < a; 

\ y(0) = 0; y(o) = 0. ^ i4j 

Similarly, the zeros of y'(a,k) exactly corresponds to the eigenvalues of the Sturm-Liouville problem 

' y" + (k 2 n{r) - $)y = 0, < r < a; . . 

y(0)=0;y'(a) = 0. 



Under the transform (3.11 1, the zeros of z(N(a),k) exactly corresponds to the eigenvalues of the Sturm- 
Liouville problem 

/ z" + [k 2 -p(O]z = 0, 0<£<N(a); 

\ z(0) = 0; z(N(a)) = 0. [6 ° j 

Similarly, the zeros of z'(N(a), k) exactly corresponds to the eigenvalues of the Sturm-Liouville problem 

z" + [k 2 -p(0]z = 0, 0<£<N(a); 

z(0) = 0; z'(N(a)) = 0. [6 ' > 



Under transformation (3.11), the quantity z 1 (N(a),k) and (z l )'(N(a),k) corresponding to refraction 
index n l (r) has common zeros. Hence, the Sturm-Liouville problems 



{z 1 )" + [k 2 -p % {C)]z l = 0, < f < N(a); 
z*(0) = 0; z l {N{a)) = 0, 



have the same eigenvalues for both i = 1,2. Similarly, 



(z l Y' + {k 2 ~p l (£)}z l = 0, < £ < N(a); 
z*(0) = 0; (z l )'(7V(a)) = 0. 



(3.18) 



(3.19) 



have the same eigenvalues for both i = 1,2. Given p l g C 2 , i = 1,2, Theorem 3.1 says these imply 
p 1 = p 2 . Therefore, we have 

^(^,fc) = 2 2 (e,fc),Ve,fc. (3.20) 



This says that 



n 1 (r)] 4 y 1 (r, k) = [n 2 (r)] 4 y 2 (r, fc) , Vr, k 



(3.21) 



Since n l never vanishes, the solutions y 1 ( r, k), y 2 (r, fc) have common zero set in C of its maximal dens ity as 

described by Cartwright's theory. From (2.27), their density is N ^ which can be computed using (2.41 ) 
and (12.42 ). Hence, 



n 1 (ppdp= / [n 2 (p)}^dp, Vr. 



(3.22) 



We have n 1 = n 2 by assumption (1.1 ). 
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